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Abstract 

~^ ■ In this paper, we describe a general method for computing Selberg-hke integrals 

based on a formula, due to Kaneko, for Selberg-Jack integrals. The general principle 

^ I consists in expanding the integrand w.r.t. the Jack basis, which is obtained by a 

Gram-Schmidt orthogonalization process. 

The resulting algorithm is not very efficient because of this decomposition. But 
for special cases, the coefficients admit a closed form. As an example, we study the 
case of the power-sums since for which the coefficients are obtained by manipulating 
generating series by means of transformations of alphabets. Furthermore, we prove 
that the integral is a rational function in the number of variables which allows us 
to study asymptotics. As an application, we investigate the asymptotic behavior 

2 ' when the integrand involves Jack polynomials and power sums. 
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^ ; 1 Introduction 

^ ! This paper is devoted to Selberg-like integrals, that is: 
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S^ . where / is a multivariate polynomial. 



For a review on Selberg integrals, see [1]. 

This integral appears in the physical problem of quantum transport. A link has been 
found between the properties of some large random matrices and the fluctuations of the 
electronic conductance in certain disordered systems (see [I], [5], [I3] or [9]). In the case of 
a system constituted by a phase coherent disordered region connected to two reservoirs of 
electrons by two ideal leads, it is possible to approximate the scattering matrix (the matrix 
that links the wave functions of the incoming and outgoing electrons) by a random matrix 
S which belongs to one of Dyson's circular ensembles. As usual, the eigenvalues of this 
random matrix are of great interest, and some experimentaly measurable variables can be 
computed with them. The contraints of the physical problem restrain the admissible form 
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for the joint probability density of the eigenvalues of S and this leads to the computation 
of integrals of the form (/)^f,K- The properties of large physical systems are related to 
the limit N ^ oo. 

In [TT], Luque and Vivo described a method for computing the asymptotics when A^ — > 
oo of {x^)ab K- Their method involves different basis of symmetric functions. Surprisingly, 
the result has an interesting combinatorial interpretation in terms of Dyck paths. 
The main tool of the algorithm presented in [11] is the expansion of power sums 

k 
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Pk{X) ■=x\ + --- + x 



in the Jack basis P^ • The coefficients arising in the decomposition 

Pk = Y,ax,kPP (2) 

Ahfc 

do not depend on the number of variables A^. Now, by the well-known formula due to 
Kaneko (see [7], 0,112]): 

ipi^yN i-r r(A, - A, + k{j -i + l)) ^ r(A. + g + K(iV - r))T{h + k{N - i)) 

^A ;aA. n r(A.-A, + /.(j-z)) \_\ T{K + a + h + K{2N-^-l)) ' 

(3) 

Luque and Vivo obtained the exact value of °' '^ with its asymptotic behavior. But 

V-)a,b,K 

they did not explain how to obtain the coefficients a\^k- 

In general. Jack polynomials are obtained by orthogonalizing the Schur basis w.r.t. a 
deformation of the usual scalar product. But this does not provide an efficient algorithm. 
We show how to improve it by computing directly the value of the coefficient a\^k- 

This paper is organised as follows. In the next section, we recall some basic properties 
of Jack and Macdonald polynomials. Section [3] contains the computation of ([T]) for Jack 
polynomials. In the last section, we give the exact value of the Selberg - power sum 
integral and present two asymptotic relations suggested by numerical results. 



2 Jack and Macdonald polynomials 

In this section, we recall some definitions and properties of symmetric functions and Jack 
and Macdonald polynomials that will be used essentially in section 14. 2[ 

2.1 Symmetric functions, Cauchy Kernel and A-ring operations 

It is well known (see e.g. [10]) that the algebra of symmetric functions is endowed with 
a A-ring structure. As usual, we will denote by (^^(X) the Cauchy function, that is the 
generating function of the complete functions Si over the alphabet X = {xi, . . . , x^}'- 

a^{X) = Y,S^{X)z^ = \{-^. (4) 



Remark that o.^ and the exponential generating series of the power sums are related as 
follows: 



a.(X)=exp('5^^.") 

\n>\ ) 



(5) 



The sum of two alphabets X and Y is defined by 

a,(X + Y) = a,(X)a,(Y) = J^ ^,(X + Y)z\ (6) 

i 

The multiplication of an alphabet by a constant is defined by extending to any complex 
number u the following property: (T2(2X) = cr^(X) which is obtained with X = Y in ([H]). 
Thus, 

a,(wX) = a,"(X). (7) 

In particular, \i u = —1, a-^—^ = ^^(X)^^. Hence, multiplicities could appear in the 

alphabet: the notion of alphabet is not restricted to sets of variables but is, by this way, 

extended to series. For example, we will use the alphabet j^, which will be considered 

as the difference of the alphabets l + t + --- + t"' + ... and u + tu + ■ ■ ■ + t^u + . . . . In 

terms of operations on symmetric functions, this is equivalent to the operation that maps 

1 -m" 
Pn on . 

^ 1 - 1" 

Finally, the product of two alphabets is defined on the complete functions: 



ai(XY) = 5^5.(XY) = nnr^ 



xyt 

The function (Ti(XY), denoted by K{X,Y), is the Cauchy kernel. It is the reproducing 
kernel associated to the usual scalar product on symmetric functions. Therefore, it satis- 
fies the following property: for two basis A\ and B\ in duality, and two alphabets X and 
Y, one has 

K(X,Y) = ^Aa(X)Sa(Y). (9) 

A 

In particular, for Schur functions Sx = det {Sx.^i+j)^ ., one has 



A'(X,Y) = 5^5a(X)Sa(Y). 



2.2 Jack polynomials 

(-) 
Jack polynomials P^ "" are defined as the unique family of symmetric functions orthogonal 

w.r.t. a one-parameter deformation (,)i of the usual scalar product such that, for two 

partitions A and /x, 

/iy(A) 

{Px,Pf,)i = zxi-j Sx,^^, (10) 

e{x) 

(where zx = Y]i"^'^^'mi{X)\ and p^ = TTp^J and such that 
i>i i=i 

pi^^(X) = mx{X) + Y,ux,,^rn^i^). (11) 

^J■<x 



where nix denote the monomial symmetric functions. 

(-) 
By definition, it is possible to construct the basis P^ by applying Gram-Schmidt 

algorithm to the Schur basis w.r.t. (, )i and beginning with the partition 1". For example, 
let us construct the py~) basis of symmetric functions of degree 3. We begin with 

r>\~) c 

-* 111 ~ '-'111. 
Then, 

J^2l — 'J21 /j_\ /j^N -Till • 

/pu; pUJ\ 

V-' 111 1 -'111 /i 

(-) (-) (-) 

The scalar products (>S'2i,iiii )i and (-Pni ) -^111 )i ^^^ known from the previous steps 
and the decomposition of Schur functions in the power sums basis. The orthogonality 
property flTUl) of the power sums gives the coefficients. 
Finally, 

3 ~ "^3 /j_x /j_x J^21 (i\ (i\ -'^111 • 

/pu; pUJv /pu; pUA 

\-'21 ) -' 21 /i V-'^lll ) -' 111 /i 

K K 

Note 2.1 Remark that the normalization property [Tl]) is verified since 

sx = mx + '^ wx,f,m^. (12) 

2.3 Macdonald polynomials 

The family of symmetric and homogeneous Macdonald Px{q,t) polynomials is uniquely 
defined by its orthogonality w.r.t. the {q, t)-deformation (, )q^t of the usual scalar product 
on symmetric functions such that 



1 - g^" 



{P\,Pf^)q,t = ^A n ^_^;,/ a,M, (13) 

j=l 
and by its dominance property: the expansion of Px{q,t) in terms of mx satisfies 

PA(X;g,t) = mxiX) + ^ t;A,^m^ (X) . (14) 

M<A 

Denoting by Kg^tO^, Y) the reproducing kernel associated to (, )q^t, one has 

i^,,i(X,Y)=ai(XY|^), (15) 

and 

K,,i(X, Y) = J] P,(X, q, t)Qx{Y, q, t), (16) 

X 
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where Qx{q, t) is the dual basis of the Px{q, t) for the scalar product (, )q^f Equivalently, 

{Pxiq,t),Q,{q,t))g,t = 6x,^. (17) 

Note that Q\ and P\ are proportional (by definition). 

Macdonald polynomials are a generalization of Jack polynomials: substituting q by t~ 



and then taking the limit t — ;■ 1, we map P\{q, t) onto P; 



(^) 



3 Exact computation of the Selberg-Jack integral 

In this section, we are interested in the computation of 

/ p{~)\N 



\ X I a,b,K. 



N 



\ I a,b,K 



18) 



It is clear from ([3j) that the asymptotics can not be precised directly since the number of 
factors of the products depends on A^. Therefore, we have to simplify this expression. 



Lemma 3.1 For any k E C, 

'i(X) iiX) 



{prnl 



b,K,N 



nn 



r(A, - A,- + Kjj - 2 + 1))T{k{j -i + 1)) 



nn 

i=l j=0 



k{N + 1-i) +j 



nn 

i=l j=0 



a + k{N -i)+j 
a + b + k{2N -i- 1) +j 



(19) 



Note 3.2 It is easy to check that we recover equation (4) of JB] by setting k = 1 in the 
previous formula. 



Proof: We need the following property of the Gamma function: 



n-l 



r{z + n) = Y[{z + i)r{z). 



(20) 



1=0 



^(A) ^(A) ^(A) AT 

First, we split the first product of (|3]) into several parts: TT TT , TT TT and 

i=i j=i+i i=i j=e{x)+i 

N N 

TT TT . The first part appears directly in (TT9|) . The last part does not depend 

i=i(X)+lj=i+l 

on the partition and gives 1 when we divide {Px /abn ^J {^)abK- 
Assuming that k G N, we can write the second part as 



nn n 

t=o j=i j=e{x)+i 



\ + i^U -i) +t 



(21) 



Let us define VT := TT TT TT ^^ + '^U ^)+\ q^^ ^^^ 
t=o i=i j=e{x)+i ^-^ ' 

To prove the relation, we begin by shifting the variable t in the numerator: t' = t + Xf 

Xi+K-l N 

«„ n n (*+-o-')) 

11 K-1 JV 

^=^ n n (i+<j-^)) 

t=o j=eix)+i 
If Aj > K — 1, we can multiply numerator and denominator by the factors for k < t < Xf. 



iiX) N 



Xi+K-l >,._! 

n iM3-i) + t)il«3-t) + t). 



w=n n -^1 s (23) 



i=i j=e{x)+i 

t=0 t=K 



' n('^(^'-^)+^) n ('^(•^' - ^) + ^) 



If Aj < K — 1, we can also decompose W following the previous equation. In both cases, 
we have 

,« n n ('+««-')) 

w = TT "" '-""'^' 

11 A,-l N 

'=' n n it+<j-^)) 

t=o j=e{x)+i 

Ai~l N 

„„n n («+''0-'+i)) 

n t=o j=e{x)+i 
Ai-l AT 

^" n n it+<j-^)) 

t=o j=e{x)+i 

using another shift of the variable t: t = t' + k. A last shift over j gives the following 
relation: 

Ai-l N+1 

^-UB^^ • (24) 

t=o j=e(x)+i 

The terms that are different in the numerator and the denominator are obtained for 
j = i{X) + 1 and A^ + 1. After simplification, we obtain (!22|) . 
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To prove (TT9l) . it is enough to use (120|) and the definition of {Px )abKN- Now, since both 
left- and right- hand sides of fl2^ are entire functions of exponential type of A^, Carlson's 
theorem |2] allows us to extend the result for all «; G C D 

From (1191) and fl2ip . it is possible to give an approximation of the asymptotic behavior 

of {Px YabKN "when N ^ oo. Indeed, {Px )ab kN ^^ ^ product of rational fractions in 
A^. Since A^ appears in the numerator and the denominator of the last product and not 
at all in the first, the asymptotic behavior depends only on the second product. In this 
product, A^ appears to the power |A|. Therefore, 

Proposition 3.3 

{Pi'hlb,.,N- N\'\. (25) 

4 Asymptotics of the Selberg - power sum integral 

4.1 General method 

/f\N 

In this section, we give the algorithm that allows the computation of {f/ahnN ~ \n ''^ 
for any polynomial /, not necessarily symmetric. First, we symmetrize /: 

®^ = ]^ I] ^-^ = ]^ S] /(a;a(i), • • • , Xa(N)), (26) 

and we remark that {&f)abKN — {f/abnN- The algorithm to compute the integral will 
then be as follows: 

(-) 

1. Expand (5/ in terms of Jack Px polynomials; 

2. Replace each occurrence of -P^ by {Px )abKN which can be computed with flT9l) . 

If (3/ admits a decomposition in the Jack basis which does not depend on the number 

of variables, the asymptotics can be exactly computed since the resulting integral is a 

rational function in A^. 

We are interested in the value of {pku bnN- ^^ ^^i^ case, we do not need to symmetrize 

(-) 
the integrand and it is sufficient to find the value of the coefficient ax,k of /\ in pk. 

This is the aim of the next section and we will find it by manipulating a specialization 

of Macdonald polynomials. 

4.2 About a specialization of Macdonald polynomials 

The aim of this section is to show how to compute the coefficient of the Macdonald 
polynomial P^' in p^. 

Proposition 4.1 With the following notations: 

• ax{s) denotes the arm-length of s (if s = {i,j), a,x{s) = Aj — j); 



• lx{s) the leg-length of s (if s = {i,j), lx{s) = X'j - i), 
the coefficient of Px{q,t) in pk is given by 

n (t'-'-^n 

(i,i)eA 
^l_qk^^;£f^ . (27) 

Proof : The coefficient of -Pa(X; q, t) in Pfc(X) equals the coefficient of c^' (X) in Q\(X; q, t) 
where c^'* = 2;^(g, t)~^j9^(X) denotes the dual basis of (p^) for the scalar product ( , )q^t 
and with 

1-q 



^(^) . -A. 



(^'^) = ^^nT3^- (28) 



i=l 



If /3^' denotes the coefficient of p^ in Q\(^; q,t), one has, because of the interpretation 
of the alphabet j^ in terms of operations on symmetric functions (see l2.ip . 



1_„ . ^^.K!-!."' 



0.(^;.,<) = E<nT^- (29) 



4=1 

Dividing left- and right- hand sides of this equation by , we can write that 

I — u 

Y3^yA(Y375?'^)-Z^/^A;^(i_tm)(i-^)llT3^- 

/i<A ^ ^^ ^4=2 

1 — m'^1 

The only term that does not vanish for u = 1 is — r because of the pole at 

^ (1 -t'^i)(l-t/) ^ 

this point. Therefore, if we take the limit u ^ 1, the partitions that give a non zero 
contribution have one part: fi = \X\ = k. Since 

1 — u^ 
lim = k, 

u^l 1 — U 

the following relation holds: 



But Qxijzf'^Qji) is known and can be expressed with help of formulae (8.3) p352 and 
(8.8) p 354 of [H]: 

n it'-' - 1'-'^ 



u 



Qx{- — -;q,t) 



1-t ' TT(1 _ ^a^{s)+mx{s)y 

sex 



JJ(1 _ qa,{s)+l^l. 
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Therefore, 



n i^'-'-i' 



sex 
and the coefficient of Pa(X; q, t) in pki^) is equal to 






sGA 



n i^-'-i'-') 



(31) 



{»,i)eA 
= (1 _ fc) (l^)f(ii^ □. 

J](l - g"-W+V^W) 

4.3 Exact computation 

The result of the previous section involves Macdonald polynomials. We have to come 
back to Jack polynomials. Applying the method described in section 12.31 (set q = t~ and 
take the limit t — )■ 1), one recovers a result of [B]: 

(-) 
Corollary 4.2 The coefficient of P^^'" in pj. is equal to 

n {{J-l)-^{^-l)) 

(»J)eA 
Oi\k = k — z=z . [61) 

\\{ax{s) + l + lx{s)K) 
sex 

Applying the algorithm, we have the following proposition: 
Proposition 4.3 

n {{J-l)-^{^-l)) 

\Pkl a,h,K,N — '^ 2-^ TT / / \ y , \ \ ^^ /a,b,K,N- I'^'^J 

xhk [[{ax{s) + l + lx{s)n) 
sex 

We have written the integral {Pk/ahnN ^^ ^ rational fraction in A^. Now, the study 
of the asymptotics is easy whence simplifying this fraction. Numerical evidences suggest 

1 a 

that, if a, h and k are constant, {pk}a,b,K '■= hm -T;{Pk)abKN ^'^^^ iiot depend on a, h 
and K and its value involves central binomial coefficients: 

{Pk}a,b,. = ^{^\. (34) 



Even when a and b are linear in A^, the hmit of -^{PkiahnN seems to converge. In 
particular, when a = k{£ — 1)N and fe = 0, we find experimentally that its value is related 
to the number of symmetric Dyck paths counted by number of peaks: 

1/ \tt f '1^'^ fk -l\fk- 1\ ^, 

In [5], we have completely described, in a combinatorial way, what happens when 
K = 1 and, in a forthcoming paper, we will prove the result for formal parameter k using 
analytic tools. 
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